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ABSTRACT 

Small devices were attached near the tip of a hov- 
ering rotor blade in order to alter the structure and 
trajectory of the trailing vortex. Stereo particle im- 
age velocimetry (PIV) images were used to quantify 
the wake behind the rotor blade during the first rev- 
olution. A procedure for analyzing the 3D-velocity 
field is presented that includes a method for account- 
ing for vortex wander. The results show that a vortex 
generator can alter the trajectory of the trailing vor- 
tex and that a major change in the size and intensity 
of the trailing vortex can be achieved by introducing 
a high level of turbulence into the core of the vortex. 

NOMENCLATURE 


a 

speed of sound, ms” 1 

A 

area enclosed by circulation path, m 2 

c 

chord length, m 

c t 

thrust coefficient, thrust/p7rft 2 R 4 

M 

hover tip Mach number, Rfia” 1 

N 

number of points to be averaged 

r 

radial distance from vortex center, mm 

r c 

vortex core radius, mm 

R 

radius of rotor tip path, m 

Re 

hover tip Reynolds number, cRQi/~ l 

U, V , w 

velocity compontents, ms -1 

Vc 

vortex convection velocity, ms” 1 

Vi 

vortex interior velocity, ms -1 

v$ 

circumferential velocity in x-y plane, ms 

VrB 

vq relative to convection velocity, ms” 1 


coordinates relative to rotor tip, mm 


mean values of x, y, mm 


x , y relative to basic blade vortex, mm 

a 

angle of position vector, deg 

0 

angle of u -f v vector, deg 


T circulation, m 2 s -1 

circulation at particular value of r, m 2 s“ 1 
S radial departure from mean location, mm 

v kinematic viscosity, m 2 s _1 

( 7 standard deviation of vortex wander, mm 

^ azimuthal angle measured from blade, deg 

w z vorticity normal to x-y plane, s” 1 

Q rotational speed of rotor, RPM 

ABBREVIATIONS 


BERP 

British Experimental Rotor Program 

BVI 

Blade-vortex interaction 

CCD 

Charged-coupled device 

FFT 

Fast Fourier transform 

LDV 

Laser Doppler velocimetry 

PIV 

Particle image velocimetry 

TG 

Turbulence generator 

VG 

Vortex generator 


INTRODUCTION 

Although the helicopter performs well in hover 
and low speed forward flight, it lacks community ac- 
ceptance in populated areas due to the noise that it 
produces. A major source of helicopter noise comes 
from the rotor blade as it cuts through its own wake. 
This phenomenon is known as blade/vortex interac- 
tion (BVI) noise, and it occurs primarily when the 
helicopter is descending. Applying a mitigating de- 
vice that reduces rotor noise can be quite challenging 
since the offending trailing vortices are shed from the 
tips of the rotor blades where the centripetal acceler- 
ation is highest. In addition to the force that such a 
device must withstand, consideration must be given 
to rotor performance since any alterations to the tip 
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region, where the dynamic pressure is highest, could 
degrade the aerodynamic efficiency of the blade. In- 
creasing the chord length and reducing the tip speed 
of the rotor may seem like an obvious solution for 
reducing BVI noise, however the added weight and 
control loads may be offsetting. 

Earlier studies have identified two important pa- 
rameters governing BVI noise generation. The first 
parameter is the miss distance between the rotor 
blade and the trailing vortex. The second parame- 
ter is the size (or intensity) of the vortex. The miss 
distance and the vortex core size can be simultane- 
ously increased with a nonplanar tip configuration 
(refs. 1,2), however, the shape of this blade limits 
its application to hover and low speed flight. Another 
approach to lowering BVI noise involves a momentary 
change in blade pitch at just the right azimuth using 
higher harmonic control (ref. 3), an active flap (ref. 
4), or individual blade control (ref. 5). Reductions 
in BVI noise have been achieved, but it is not certain 
whether these improvements are the result of an in- 
crease in the miss distance, a decrease in the intensity 
of the interacting vortex segment, or a combination 
of both effects. There is also some concern that the 
inputs needed for noise reduction may promote an 
increase in vibration (ref. 5). 

There have been numerous attempts to simply 
reduce the intensity of the trailing vortex. An early 
effort was aimed at reshaping the planform of the 
blade near the tip such that the local flow swirling 
around the streamwise edge of the blade would not 
encounter a surface on the upper side on which to 
reattach (or stagnate). This design was known as 
the ogee tip and it did result in a significant reduc- 
tion in the concentration of vorticity in the wake of 
a model wing (ref. 6). The same design was applied 
to the blades of a model rotor and was found to pro 
duce a similarly difFused tip vortex, however, there 
was an unacceptable degregation in performance at 
the higher thrust levels due to an early onset of stall 
(ref. 7). Other attempts, such as tapering the plan- 
form or the thickness, sweeping the tip, or adopting a 
major planform change (BERP), have all resulted in 
only modest increases in the diffusion of the vortex, 
with the conclusion that more intrusive devices need 
to be considered (ref. 8). 

Adding devices to the tip of the rotor is represen- 
tative of the more agressive measures taken to reduce 
BVI noise. The first known attempt originated as a 
spoiler placed on the upper surface near the wing tip 
of a full-scale transport (ref. 9). Although the maxi- 
mum swirl velocity was reduced by a factor of 3, there 
was little difference in the rolling moment experienced 
by the aircraft following in the wake of the trans- 
port (ref. 10). Nevertheless, the dramatic increase in 


core size was sufficient to attract its application to a 
model rotor (ref. 11). The device, which functioned 
as a spoiler, was placed at the quarter chord, nor- 
mal to the upper surface, and measured 9% in height 
and 12% in width relative to the chord of the rotor 
blade. Tests revealed a substantial reduction in BVI 
noise, however, a discouragingly high level of power 
was consumed and there was a significant increase in 
broadband noise (ref. 11). 

One of the less intrusive additions to the tip of 
the rotor blade was the subwing (ref. 12). This de- 
vice took the form of a small wing with a chord length 
about 20% that of the rotor blade. The subwing be- 
came an extension of the rotor blade and was intended 
to force the circulation near the tip of the blade to 
be shed as two co-rotating vortices, each having less 
intensity than would occur without the subwing. Re- 
sults from this test showed a slight reduction in torque 
at moderate thrust levels and a modest reduction in 
the swirl velocity of the vortex trailing from the ro- 
tor blade. A more recent test revealed that the two 
vortices initially contract more rapidly than for the 
blade-alone vortex, but then combine into a single 
vortex (at around ip = 150°, depending on the lift on 
the subwing) and follow the same trajectory observed 
for the blade-along vortex (ref. 13). 

Recognizing the effectiveness of a spoiler in dif- 
fusing the trailing vortex, various devices have been 
attached to the tip of a fixed wing at locations both 
near and beyond the trailing edge (ref. 14). Once 
again, several designs proved quite capable of diffus- 
ing the trailing vortex and each carried a certain mea- 
sure of drag penalty. Since any device that causes an 
increase in drag on a wing would be expected to be 
even more objectionable when applied to the tip of 
a rotor blade, the less intrusive among these spoilers 
have been tested on a model rotor (ref. 15). Small an- 
gles and wedges were attached to the trailing edge of 
the rotor blade, and a sizeable reduction in BVI noise 
was observed. Although the results were encourag- 
ing, the seemingly inevitable increase in power led to 
the suggestion that some means of actively deploying 
these devices might be the only reasonable solution 
(refs. 8, 11, 15, 16). 

While avoidance, rather than forced diffusion, 
may be the best approach for fixed-wing aircraft to 
deal with the hazards of across-trail and in-trail vor- 
tex encounters (ref. 10), rotor blades engage rela- 
tively young vortices so that methods that even tem- 
porarily enlarge the trailing vortex could prove useful. 
Although the concepts of vortex pairing and turbu- 
lent diffusion are not new, the present study examines 
these ideas on a rotor and with a diagnostic technique 
that provides a level of detail that was not previously 
available. 
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Similar to the wake of a fixed-wing aircraft, the 
dominant feature of a rotor wake is the trailing vortex 
that is shed from the tip of the rotor blade. Experi- 
ments typically focus on the location, size, and inten- 
sity (or magnitude of vorticity) of the trailing vortex 
at various wake ages. The trajectory of the vortex 
can be accurately determined from flow visualization 
based on a triangulation technique using laser light 
sheets (ref. 17). While this technique is relatively 
easy to apply, it does not provide any information 
about the structure of the vortex. Although tedious, 
the structure of the trailing vortex can be obtained 
using laser Doppler velocimetry (LDV). It is impor- 
tant to note that small disturbances in the freestream 
flow will cause the vortex to wander, or meander, over 
a lateral distance that increases with distance from 
the blade. A point measurement technique, such as 
LDV, will therefore require some procedure for con- 
ditionally averaging the results (ref. 18), especially 
as the wake ages or when the freestream turbulence 
level is high. A planar measurement technique, such 
as PIV, has the advantage of requiring substantially 
less time to acquire the 3D velocity field. Neverthe- 
less, this data must also be conditionally averaged 
when vortex wander is present (ref. 19). 

This paper contains a discussion of the test setup 
and a method for obtaining the 3D velocity field in 
the wake of the rotor at specific wake ages. The proce- 
dure for identifying separately the centers of vorticity 
and swirl, and how they relate respectively to condi- 
tional ensemble averaging and the vortex structure, 
will be presented. Data will be shown that clearly 
expose the distinguishing features of the flow and the 
effectiveness of both vortex pairing (with same and 
opposite sense) and forced turbulent diffusion on the 
aging trailing vortex. 

TEST DESCRIPTION 

Test Chamber and Stand- The experiment 
was performed in the Hover Test Chamber at the 
NASA Ames Research Center under the authority of 
the U.S. Army Aeroflightdynamics Directorate. The 
chamber has a base of 26 ft x 32 ft and a height of 
28 ft. To minimize recirculation in the chamber and 
to limit the influence of the floor on the flow, the col- 
lective pitch of the rotor blade was set to a negative 
angle so that the wake would be directed upwards. 
With the rotor acting like a pump, air was drawn into 
the chamber through filters placed across two oppos- 
ing rollup doors, then confined by an annular diffuser 
located above the rotor, and finally exhausted to the 
exterior through openings near the top of the chamber 

(%• I)- 

The rotary-wing test stand that was used to drive 
the rotor was configured with a single, 90-hp electric 


motor and a 2.5:1 transmission. A flexible coupling 
between the input shaft and the dummy balance was 
instrumented to measure torque. An encoder with a 
resolution of 4096 steps per revolution was attached 
to the rotor shaft. The encoder signals were passed to 
a variable delay circuit that enabled the PIV cameras 
and Nd:YAG laser to be synchronized to any desired 
rotor azimuth. 

Rotor Blade and Hub- The rotor consists of 
a single counter-weighted aluminum blade having a 
rectangular planform, zero twist, 7.5 in. chord, and 
a 45 in. radius (fig. 2). The outer 50 percent of the 
blade conforms to a NACA 0012 profile (thickness at 
the tip is 22.9 mm). Over the inner portion of the 
blade radius the profile linearly thickens to a NACA 
0020. The rotor hub consists of two steel sections 
that are clamped together to hold the blade at a 0° 
coning angle. The collective pitch angle was fixed 
at -8°. When mounted on the test stand, the rotor 
and hub were approximately 7.5 ft above the floor of 
the test chamber. The solidity of this single-blade 
configuration is 0.053. 

Blade Attachments- One of the mechanical de- 
vices under study was designed to function as a vor- 
tex generator and the other as a turbulence generator 
(figs. 2 and 3). The vortex generator has a NACA 
0012 profile, rectangular planform with a rounded tip, 
and measures 0.12c in chord and 0.10c in span (where 
c refers to the chord length of the main element). 
The quarter-chord axis of the vortex generator passes 
through the quarter-chord location of the main ele- 
ment. The turbulence generator consists of a flat rect- 
angular section and measures 0.12c in length (with a 
fixed alignment with the span of the main element) 
and 0.10c in height (oriented normal to the surface 
of the main element). The center of the turbulence 
generator also passed through the quarter-chord axis 
of the main element. Both generators were mounted 
0.057c inboard from the tip of the main element. If 
the vortex generator were to be rotated 90°, its pro- 
jected area would appear to the oncoming flow the 
same as the turbulence generator. The relative di- 
mensions and placement of the turbulence generator 
are similar to the spoiler employed in an earlier study 
(ref. 11). 

Test Conditions- The present test was per- 
formed at a constant rotor speed of 870 rpm (14.5 
Hz), which corresponds to a tip speed of 341.7 ft/sec. 
Based on an average ambient temperature of 65° F 
and barometric pressure of 760 mm Hg (14.7 psi), the 
Reynolds number (based on c) at the tip of the ro- 
tor blade was R e — 1.33 x 10 6 and the Mach number 
was M — 0.30. Based on a collective pitch angle of 
-8°, the equivalent two-blade thrust coefficient was 
Ct = 0.005 (ref. 20). The vortex generator was tested 
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at eight angle-of- attack (relative to the chord of the 
main element) conditions: -15°, —10°, — 5°, 0°, 5°, 
10° , 15°, and unrestrained. The turbulence generator 
was tested at a fixed orientation normal to the up- 
per surface of the main element and parallel to the 
leading edge of the blade. 

Stereoscopic PIV Concept- The velocity of a 
fluid is inferred from the motion of discrete particles 
that are suspended in the flow. A thin sheet of light is 
used to define the plane of interest. The particle im- 
ages that are recorded at two different times are then 
cross correlated to yield displacement vectors (mag- 
nitude and direction). The time interval between the 
two images must be short enough that the particles 
remain in the light sheet, yet long enough that the 
particle displacements are perceptible. It is assumed 
that the particles are small enough to accurately track 
the flow. 

A single camera placed normal to the sheet of 
light results in a two-dimensional array of velocity 
vectors. If the flow has a significant third com- 
ponent (directed normal to the plane of the light 
sheet), an error due to perspective develops in the 
two-dimensional array that is zero in the center and 
increases toward the image boundaries. This seem- 
ingly undesirable sensitivity to out-of-plane motion 
was later exploited to derive the third component of 
velocity (ref. 21). In fact, sensitivity to the out-of- 
plane component of velocity over the entire image 
plane can be increased by intentionally placing the 
camera at an oblique angle to the light sheet. To 
uniquely determine a three dimensional particle dis- 
placement requires two cameras that are oblique to 
the light sheet, each offering a different perspective of 
the particle motion between the two exposures (fig. 
4). This procedure is known as Stereoscopic PIV or 
3D-PIV. 

The technique demands both critical focus and 
maximum brightness of the particle image. Conven- 
tional lens mounts create a plane of focus in object 
space that is parallel with both the image plane and 
the lens plane. When the lens plane is oblique to 
the object plane (as is the case for Stereoscopic PIV), 
particles will be increasingly defocused as their dis- 
tance from the centerline of the object plane increases. 
The conventional corrective action would be to in- 
crease the depth of field by reducing the aperture of 
the lens. Unfortunately, reducing the lens aperture 
drastically reduces both the particle image intensity 
and the image resolution. The problem is solved with 
Scheimpflug focusing (refs. 22 and 23). To keep the 
object plane in focus when the lens plane is rotated 
to an oblique orientation, the image plane (where the 
camera sensor is placed) must also be rotated such 
that all three planes (object, lens, and image) inter- 


sect along a common line (fig. 5). 

Arrangement of Equipment- The cameras 
were mounted in a forward-scatter position on a hor- 
izontal plane that passed 3 in. above the rotor disc 
(fig. 6). The bisecting angle between the light sheet 
and each camera was 36°, and the cameras were ro- 
tated so that the image areas were coincident. The 
centerline dimensions of the image area were 14 in. 
vertical and 18 in. horizontal. Due to perspective the 
vertical dimension of the image ranged from 12 in. 
(near side to the cameras) to 16 in. (far side). The 
lower, outboard corner of the image area was placed 
so as to capture the trailing vortex over a maximum 
range of wake ages (which was less than one revolu- 
tion for this test). 

Camera Specifications- Images were acquired 
with 8-bit charged-coupled device (CCD), cross- 
correlation cameras having a sensor array of 1008 
(horizontal) x 1018 (vertical) pixels. Each pixel 
measures 9 pm on a side. The cameras can oper- 
ate in a double-exposure mode and acquire two non- 
interlaced, full-frame images in a single frame interval. 
The time interval between images is variable between 
2 ps and 30 ms. These are non standard video cam- 
eras that can be externally triggered and driven at 
any frequency up to 15 Hz in double-exposure mode. 
A computer interface provides control over gain, con- 
trast, black level, and trigger mode. Both cameras 
used 55-mm fl.2 lenses. The lenses were remotely 
translated to focus on the centerline of the image area 
and the sensor (located inside the camera body) was 
rotated about its centerline to satisfy the Scheimpflug 
condition. 

Laser and Sheet Optics- A Nd:YAG laser was 
frequency doubled to provide a beam having a wave- 
length of 532 nm, a pulse width of 9 ns, and power of 
350 mJ. The laser model used has a beam diameter 
of 9 mm and a divergence of 0.50 mrad. All optical 
elements have damage thresholds exceeding 1 J/cm 2 . 
A laser-light sheet was formed using lenses that thin 
the beam in one direction and expand the beam in 
the orthogonal direction. The thickness (horizontal 
direction in this test) of the light sheet was controlled 
by two cylindrical lenses, the focal length of one was 
f — H-200 mm and the other f — -200 mm. By adjust- 
ing the distance between these two lenses, the beam 
waist 1 mm thick) could be positioned in the re- 
gion imaged by the cameras. The vertical expansion 
of the light sheet was controlled by a single f = -75 
mm cylindrical lens. The distance between the sheet 
forming optics and the imaged area was about 25 ft. 

System Alignment- The rotor blade was posi- 
tioned at the 0° reference azimuth and a calibrated 
laser level (sweeping type) was used to establish the 
vertical plane passing across the trailing edge of the 
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blade. The laser-light sheet was adjusted to be coin- 
cident with the laser level. A calibration target was 
placed against the trailing edge, then translated to 
coincide with the area viewed by the cameras, and 
finally leveled to within ±0.01° with a digital incli- 
nometer. The location of the blade tip relative to 
a point on the calibration target was measured to 
within ±1/32 in. The alignment of the target was 
considered to be satisfactory when the surface of the 
target was evenly grazed by the light sheet. Using 
only white light to illuminate both sides of the tar- 
get, the cameras were focused and calibration images 
were recorded. 

PIV Software- This experiment was performed 
using the Integrated Design Tools (IDT) WinVu v5.10 
software (ref. 24). This software functions as both a 
data acquisition interface and an image-processing in- 
terface. This has the advantage of providing a fully 
integrated calibration procedure that allows for acqui- 
sition, quality assurance and data reduction. The cal- 
ibration yields all the optical parameters required for 
accurate reconstruction of three-component velocity- 
vector fields. A calibration was performed before and 
after each series of runs. 

Images are processed by first covering the region 
of interest with an interrogation grid. The intersec- 
tions of the vertical and horizontal grid lines define 
the centers of each interrogation window. The inter- 
rogation window is the smaller region in the reference 
and delayed images that are cross-correlated using an 
optimized Fourier transform. When the camera cal- 
ibrations are applied, the grid in each camera view 
covers the exact same area in the flow field. 

The size of the interrogation window is based 
on the particle density and the maximum instanta- 
neous displacement of the particles. Typically the 
window size ranges from 8 to 64 pixels on a side. 
Each interrogation window yields one vector. The 
cross-correlation of the interrogation window in the 
reference image with that of the delayed image yields 
a correlation map. The location of the correlation 
peak in space determines the local displacement in 
both magnitude and direction. 

WinVu incorporates quality-assurance tests in 
the vector calculations. The processing begins with 
a first-pass correlation between the reference image 
and delayed image to determine the maximum dis- 
placement range. For the second pass the software 
enlarges the interrogation area in the delayed image. 
The amount of the enlargement is based on the maxi- 
mum displacement range determined in the first pass. 
Since the Fourier transform requires both interroga- 
tion areas to be identical, the software “adds zeros” to 
the interrogation area of the reference image to match 
the size of the delayed image area. This technique 


maximizes the probability of correlating all the parti- 
cles found in the reference image interrogation area to 
those found in the delayed image, thereby maximizing 
the statistical accuracy. 

The software then counts the number of parti- 
cle images (each composed of several contiguous pix- 
els) in the interrogation window. Ten particle pairs 
are required to contribute to a correlation map. The 
centroid of each particle is calculated, thus yielding 
sub-pixel accuracy of their position. A second-order 
curve fit is determined from the ten displacements. 
This curve fit yields a single vector, at the precise grid 
point location, whose error is reduced by a factor of 
0.3 (1 /y/N, where N = 10 particles) over the straight 
fast Fourier transform (FFT) of the same area (ref. 
24). 

If there are not enough particles in a given inter- 
rogation region, the program will automatically en- 
large the interrogation area of the reference image. If 
there are ten or more particles in this enlarged area, 
then the correlation proceeds in the same manner de- 
scribed above. The data for this grid point is con- 
sidered to be a recalculation. If ten particles are not 
counted, the software will expand the interrogation 
area incrementally. It will repeat this process until 
the interrogation area expands to a limit of 64x64 
pixels. If the ten-particle threshold is not met for the 
largest area, then velocity values are interpolated for 
that grid point using a second-order curve fit based on 
nearest-neighbor values. This vector is also counted 
as a recalculation. The software tracks the number of 
these recalculations and displays that number after 
each camera view is processed. If the number of re- 
calculations is less than 1% of the total vectors, then 
the data are considered reliable. 

Calibration Procedure- The stereoscopic im- 
ages recorded by the cameras that are oblique to the 
object field must subsequently be corrected for both 
magnification and perspective. One method by which 
this can be accomplished results from recording the 
images of a fiat rectangular target (with precisely 
known dimensions) that is placed in the field of view 
of both cameras. The target used in this test con- 
sisted of a double sided print containing three rectan- 
gles (the largest measuring about 217 mm on a side) 
and a background dot pattern that simulated a par- 
ticle field (fig. 7). 

A reliable method for assessing the accuracy of 
the velocity measurements obtained from test data 
is to evaluate the displacements recorded for known 
target translations (ref. 25). Consecutive transla- 
tions of the target over orthogonal distances of 0.200 
in. resulted in an in-plane standard deviation of 0.3% 
(horizontal) and 0.5% (vertical), and an out-of-plane 
standard deviation of 1.2%. Since actual particle dis- 
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placements will normally be substantially smaller, the 
probable error should be increased by several times. 
Therefore, in the worst case, the out-of-plane com- 
ponent is estimated to be accurate to within about 
3%. 

Particle Seeding- Proper seeding of the flow 
is critical to accurate PIV measurements. The seed 
particles must be evenly distributed and of sufficient 
density to define the flow without altering its physical 
properties. The particles must also be small enough 
that they accurately follow the flow (especially chal- 
lenging in accelerating flows), yet large enough that 
they scatter a sufficient amount of light to be de- 
tected. It is usually convenient to have remote control 
over the delivery of the seed material into the flow. 
The particle generator used in this test employs an in- 
ert gas to atomize a non-toxic, pharmaceutical-grade 
mineral oil. The mist is vaporized and then condensed 
before being released into the flow. The particle size 
is estimated to be less than 0.5 //m. The particle gen- 
erator was located near the floor of the hover chamber 
so that the particles would mix with the air entering 
into the chamber before being drawn into the wake of 
the rotor. 

Data Acquisition- The CCD cameras were 
connected to separate frame-grabber cards installed 
in a PC workstation. Using WinVu software, images 
were acquired and immediately processed in order to 
evaluate the quality of the raw images and the ade- 
quacy of the inter-pulse time delay. Good image cor- 
relations depend on such issues as background light 
contamination, particle image brightness, contrast, 
focus, beam alignment, and light pulse separation. 

Prior to taking data, fine focusing and beam 
alignment were performed using a seeded jet. Fo- 
cusing was accomplished by translating the camera 
lenses while observing the real-time image displays. 
The gain and black levels were adjusted to maximize 
the contrast and brightness of the particle images. 
The two laser sheets were judged to be coplanar when 
both laser pulses were recorded as a single image and 
the resulting particle images appeared as doublets. 

The correct inter-pulse delay can best be deter- 
mined under actual test conditions. Images are ac- 
quired and correlated to determine the maximum par- 
ticle displacement. This displacement corresponds to 
the maximum particle velocity and the pulse sepa- 
ration should be adjusted to produce a ± 3 pixel 
range, regardless of the window size. Pixel displace- 
ments that are too small will limit the dynamic range 
whereas displacements that are too large will decrease 
the probability of correlation (ref. 26). Because of 
the high in-plane particle displacements, caused by 
the high circumferential velocity of the trailing vor- 
tex from the rotor at different wake ages, the optimum 


pulse separations varied from 30 to 50 //s (Table 1). It 
should be noted that if the thickness of the light sheet 
is 1 mm and the minimum delay between laser pulses 
is 30 /is, the highest out-of-plane velocity (u;) that can 
be detected is 33 m/s. Since the tip of the rotor blade 
is moving at 104 m/s, portions of the flow in the near 
wake may be in error. If this condition were to occur 
at a calculation node, an interpolated value based on 
values at neighboring nodes would result. Therefore, 
at any point in the flow where w reaches 33 m/s, it 
may be assumed that the correct value may actually 
be much higher. Once image quality was assured, 50 
image pairs were acquired per wake age. Data for each 
wake age required about 200 megabytes of storage. 

Table 1: Pulse Delay Times (/is) 


Wake Age 


Config 

2° 

10° 

o 

O 

CO 

60° 

100° 

150° 

210° 

280° 

Basic 

30 

30 

30 

30 

30 

30 

30 

35 

Free VG 

30 

30 

30 

30 

30 

30 

30 

30 

0° VG 

30 

30 

30 

30 

30 

35 

35 

35 

+5° VG 

30 

30 

30 

30 

30 

35 

35 

35 

-5° VG 

30 

30 

30 

30 

30 

35 

35 

35 

+10° VG 

30 

30 

30 

30 

30 

35 

35 

35 

-10° VG 

30 

30 

30 

30 

30 

35 

35 

35 

+15° VG 

30 

30 

30 

30 

30 

30 

30 

30 

-15° VG 

30 

30 

30 

30 

30 

35 

35 

35 

90° TG 

30 

30 

30 

30 

40 

45 

50 

50 


Post-Test Data Processing- After completion 
of the test, the data was processed on a PC worksta- 
tion using Provision software, which is a more ad- 
vanced release that replaces WinVu. To obtain good 
resolution of the primary vortex, a calculation mesh 
having 99 nodes (horizontal) by 91 nodes (vertical) 
was constructed over the region of significant inter- 
est in the flow field (fig. 8). The area covered by 
the mesh resulted in an average of 8.5 pixels between 
nodes in both directions. The interrogation window 
was set at 20 pixels on a side, which gave about a 
59% overlap. The total physical area measured about 
386 mm (with horizontal increment of 3.9 mm) by 237 
mm (with vertical increment of 2.6 mm). 

A single file for each wake age includes a flag for 
every velocity measurement that indicates the nature 
of that vector (such as valid, invalid, interpolated, 
recalculated, or not calculated). Reduced data files 
containing the coordinates of each calculation node 
and the three components of velocity (in terms of dis- 
placement) were stored in ASCII format (about 32 
megabytes in size) and then transferred to a main- 
frame computer for analysis. 

The velocity components were first converted 
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from displacement units to velocity units based on 
the pulse duration for that particular measurement. 
The physical coordinates were transformed so that 
x = y = 0 would correspond to the tip of the rotor 
blade (fig. 9). The sequence for extracting the char- 
acteristics of the trailing vortex (fig. 10) for each case 
begins with a calculation of the vorticity field, for 
every image pair. Vorticity can either be calculated 
by differentiation or integration according to: 


or 

u> z — A~ X T = A~ l j) (UyV) ■ d(x , y) (2) 

where the direction of integration is such that the 
enclosed area, A, is on the left of the integration path. 
In this study the vorticity was calculated using the 
integration (or circulation box) method. 

When more than one vorticity center is present in 
the field of view, a single image pair must be processed 
to determine where the field should be partitioned so 
that statistics are built on only one vortex. A sim- 
ple two-part partition is usually sufficient to separate 
the two vortices. However, when a vortex generator 
is installed on the rotor blade, a small (but poten- 
tially intense) satellite vortex develops along side of 
the main vortex. In this case it is necessary to focus 
the interrogation on a more restricted region of the 
flow field. In such cases the statement “account for 
wander: focus” will appear in the figure for this data. 

After calculating the vorticity distribution over 
the entire image, a search was performed to discover 
the zone of greatest vorticity concentration based on 
values exceeding 50% of the maximum (fig. 11 shows 
an example of the search). The center of vorticity 
for a given image was defined to be at the average 
of these selected locations. The resulting centers of 
vorticity (one for each image) formed a set for which 
the mean value (x, y) and the standard deviation (a) 
could be calculated. Any image with a vortex cen- 
ter (x, y) that deviated from the mean (x, y) by more 
than 1.5(7 was purged from the set. The 50% thresh- 
old criterion prevents the center of vorticity for an 
image from being determined by a single, and possi- 
bly erroneous, measurement. The standard deviation 
criterion prevents an atypical departure, albeit a le- 
gitimate image, from being factored into the average. 
The threshold and purging criteria both have the ef- 
fect of reducing the degree of wander in the data (fig. 
12). When the acceptance threshold is 100% (the vor- 
ticity center in each image is determined by a single 
point), 47 points are retained and 3 points are dis- 
carded because they exceeded the standard deviation 


criterion. When the acceptance threshold is 50% (the 
vorticity center in each image is determined by an ag- 
gregate of points), a distinct population of 44 points 
results (some of which are coincident). The quantity 
S that appears in this figure, as well as many others, 
represents the degree to which the vortex wanders for 
a given set of images, and is calculated by taking the 
square root of the sum of the squares of the x and y 
offsets from the mean. 

The images must be averaged to smooth out 
small irregularities in the flow. However, whenever 
there are large features in the flow, such as a trailing 
vortex, and the vortex appears in different positions 
from image to image (vortex wander), important de- 
tails of the structure tend to be smeared out if a sim- 
ple average of the images is performed. The remedy 
is to perform the average after artificially aligning the 
images based on a recurring feature of the flow, which 
in this case is the center of vorticity. This process, in 
contrast to a simple average, is referred to as condi- 
tional ensemble averaging, the condition here being 
the alignment of all the images in the set based on 
their centers of vorticity. The image with a vortex 
center closest to the mean location that was estab- 
lished for the set of images was selected as the “an- 
chor”. The indices of each image matrix were then 
adjusted according to the offset of each vortex cen- 
ter from the vortex center in the anchor image. All 
data in the anchor image were retained. However, 
some portion of all other images with adjusted indices 
that fall outside the boundaries of the anchored image 
were necessarily discarded. Hence, as a consequence 
of vortex wander, the population contributing to the 
conditional average was greatest over the interior of 
the matrix. 

For reasons that will become clear in the follow- 
ing discussion, emphasis will be placed on the swirling 
nature of a vortex flow. To extract information about 
the geometry of the vortex, the center of the vortex 
was assumed to be at the center of swirl (which is not 
necessarily at the center of vorticity that is accreted 
from the rotor wake as a spiraling sheet with varying 
vorticity). Using the center of vorticity as a start- 
ing point, the surrounding locations are interrogated 
to determine the best node for which the sum of the 
dot products of two unit vectors is a minimum over 
a neighborhood of locations surrounding the candi- 
date node. One of the unit vectors is defined by the 
coordinates of the neighboring node relative to the 
candidate node. The other unit vector is defined by 
the velocity vector at the neighboring node. Denoting 
the angles of these two unit vectors by a and fd (fig. 
13), then the node nearest the center of swirl is found 
when the following expression is a minimum: 
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N 

^|cos(a-/?)|// V (3) 

n= 2 

where N is the total number of neighboring points 
considered. Relative to the center of swirl for a pure 
vortex in a stationary flow, a and /? will be orthogonal 
for all neighboring points and the sum will be exactly 
zero. For most real flows the sum will not be precisely 
zero, and the procedure works best when the neigh- 
borhood under consideration does not extend beyond 
about one or two core diameters (which is normally 
sufficient to cover a majority of the vorticity in a sin- 
gle trailing vortex). After locating the center of swirl, 
the size of the vortex core can be estimated. This is 
accomplished by dividing the neighborhood into an- 
nular zones and then determining the average swirl 
(or circumferential) velocity for each zone: 

AT 

X>|/* ( 4 ) 

n = l 

where N is the number of nodes within a given an- 
nular zone. The mean diameter for the zone having 
the highest average swirl velocity becomes the esti- 
mated diameter for the vortex core. This procedure 
appears to work well even though the flow is not ax- 
isymmetric. The swirl velocity, v$, is defined by the 
cross product of the unit position vector (relative to 
the center of swirl) associated with a node and the 
inplane velocity vector u -f v. This relation can be 
expressed as: 


ve = |(u + t?) sin(a - f})\ (5) 

where it is understood that vg is orthogonal to r. 

The ultimate method for calculating the size of 
the vortex core requires that the space surrounding 
the center of swirl again be divided into annular zones. 
This time a larger array of bins is constructed that 
contains the average swirl velocity and the associated 
mean radius of each annulus. This array is then fit 
with a least-squares spline subject to the condition 
that the resulting curve have only one inflection along 
its inner extent (nominally set to twice the estimated 
core radius that was previously found). The radius 
of the vortex core is now defined by the point of in- 
flection since that is where a maximum value for vg is 
reached. This procedure is considered to offer a more 
rational approach for determining the size of the vor- 
tex, especially when the velocity peaks that are char- 
acteristic of vortices (and upon which core sizes are 


traditionally based) are dependent on how the vortex 
is sliced (vertical, horizontal, or otherwise). 

The final processing stage involves the construc- 
tion of contour plots for t>, w, and together 
with profiles of these variables along a horizontal cut 
(y= constant) through the center of swirl. The convec- 
tion of the vortex does not affect w and u ) z , however 
it does produce an offset (or distortion if the convec- 
tion field is not uniform) in the u and v components of 
velocity. In this study the vortex is assumed to be uni- 
formly convected. Two methods for discovering the 
convection velocity are considered. One method as- 
sumes that the convection velocity at any given wake 
age is determined by the interior velocity measured 
at the center of swirl (the y component is defined as 
Uj). In the other method the convection velocity of 
the vortex is deduced from the record of swirl center 
coordinates calculated at each wake age (the y com- 
ponent is defined as v c ). The velocity profile for v 
relative to the convection velocity as determined by 
each of these methods (either v — Vi or v — v c ) is in- 
cluded in the presentation of results. 

RESULTS AND DISCUSSION 

Including the basic rotor case, a total of 10 con- 
figurations were analyzed at 8 different wake ages (Ta- 
bles 2 and 3). 

Table 2: Configurations 


Basic 

Free Vortex Generator 
0° Vortex Generator 
+5° Vortex Generator 
—5° Vortex Generator 
-f 10° Vortex Generator 

— 10° Vortex Generator 
-(-15° Vortex Generator 

— 15° Vortex Generator 
90° Turbulence Generator 


Table 3: Wake Ages 


2 ° 

10 ° 

30° 

60° 

100 ° 

150° 

210 ° 

280° 


The vortex generator is defined to be at a positive 
angle when the trailing edge swings towards the tip 
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of the rotor blade. For positive angles, the vortex 
that trails from the tip of the generator will have the 
same sense as the primary vortex that trails from the 
tip of the rotor blade. The actual angle of attack of 
the vortex generator may be quite different from its 
geometric angle because of the influence of the flow 
around the tip of the rotor blade on the local flow 
approaching the generator. All data were obtained 
with the rotor collective angle fixed at -8°, and all 
wake data were obtained at a constant rotor speed of 
870 RPM. 

Torque Summary- The rotor torque was 
recorded at several intermediate speeds prior to reach- 
ing the ultimate test speed of 870 RPM. These mea- 
surements are summarized in fig. 14 and indicate that 
all of the vortex generators that were oriented at fixed 
angles (0° , ±5°, ±10°, ±15°) resulted in torques that 
were lower than for the basic blade case. However, 
the “free” (or unrestrained) vortex generator shows 
an 8% increase in torque while the turbulence gen- 
erator shows a 18% increase in torque. Given that 
these devices are located near the tip of the rotor 
blade, even small changes in the drag force on these 
elements should produce a noticeable change in the 
torque. The sizable increase in torque for the turbu- 
lence generator case is therefore in keeping with the 
expected level of increase in drag for such a blunt 
configuration, and demonstrates why such a device 
should be retracted when not required. If properly 
designed, the free vortex generator would be expected 
to react to the local flow much like a weather vane, 
and therefore generate a minimum level of drag. Since 
this does not appear to be the case, the center of grav- 
ity may have been offset from the pitch axis, thereby 
rendering the element at an angle sustained by a cen- 
tripetal force that increases with RPM, or the element 
may be in a state of flutter. 

Wake Analysis- For the 50 image pairs that 
were recorded at each wake age, approximately 12% 
were discarded due to excessive vortex wander from 
the mean location. Even when the vortex wanders 
by no more than one core diameter, a conditional 
ensemble average is required to obtain an accurate 
profile of the vortex. In contrast to a preferred con- 
ditional average of the data, a simple average results 
in a lower peak velocity and an increased core radius 
(fig. 15), although the circulation calculated along a 
circular path in the outer 1/r region appears to be 
independent of the method of averaging. The same 
conclusions have been reported using an exponential 
expression for the circumferential velocity and sim- 
ulating different amplitudes of vortex displacement 
(ref. 27). It was shown that a substantial degreda- 
tion in the velocity profile occurs when the amplitude 
is greater than one core radius. 


Circulation, T, carried by the entire image will 
be presented for each case. This value is obtained by 
summing all of the individual circulation-box calcu- 
lations within a given image, which amounts to 8,820 
contributions for the grid used in this analysis. The 
total image circulation is quite sensitive to the con- 
vection of vorticity through the scene. In very general 
terms, the circulation calculated for a scene that con- 
tains slices of the vortex from two revolutions would 
be expected to decrease by about 50% when the sec- 
ond slice is convected out of the scene. The total 
image circulation is useful when comparing the wakes 
from different rotors (or with different vortex diffu- 
sion devices) and for evaluating the overall accuracy 
of numerical codes. 

For an “ideal” vortex, all of the vorticity is con- 
tained within r c and the circulation remains constant 
for r > r c . However, in the case of a “real” vortex 
that is isolated and two-dimensional, the circulation 
within this body of fluid will rapidly increase with r 
until reaching some point beyond the core radius. For 
larger values of r the circulation will approach a con- 
stant value as long as the integration path remains in 
the outer “inviscid 1/r” region. For the present rotor 
wake case the circulation appears to follow the clas- 
sical trend and approach an asymptotic value only in 
the immediate vicinity of the vortex, but then devi- 
ates from a constant value at distances greater than 
about twice the core radius (fig. 16). This devia- 
tion is even less satisfactory after taking into account 
the vortex convection velocity. Since the axis of the 
trailing vortex from a rotor actually follows a heli- 
cal path, the flow field around any given segment of 
the vortex will be affected by the induced flow caused 
by neighboring segments from previous revolutions of 
the blade. Another factor to be considered is the con- 
tinual (although diminishing) accretion of chordwise 
vorticity that was initially shed into the wake along 
the span of the blade due to changing circulation. Re- 
garding the variation of vq with r, there is less scat- 
ter in the relative-velocity data than there is in the 
absolute-velocity data. This is because in the rela- 
tive case the flow field is referenced to the convection 
velocity of the vortex, which renders it more axisym- 
metric. 

Given that any experimental technique is subject 
to some degree of error, it is useful to compare the 
present results with earlier measurements obtained 
under similar conditions. After applying a coordinate 
offset and accounting for the different tip speeds, pre- 
viously obtained LV data (ref. 20) can be compared to 
the present PIV data (fig. 17). The PIV data shows 
a shortfall in the magnitude of the velocity peaks as 
well as an overestimate of the core size. Since the 
data in both cases was taken in the near wake, vor- 
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tex wander (which was not accounted for in the LV 
data) is not believed to be a factor (which if any- 
thing, would have lowered the peak velocities in the 
LV data). Assuming that the LV data is more cor- 
rect, the difference is most likely the result of a PIV 
inter-pulse delay time that was too long to capture 
the particles having a higher than anticipated out-of- 
plane component. When this happens, the reported 
velocities at the grid points where this may have oc- 
curred would have been the result of interpolations 
using lower values at neighboring grid points. As ex- 
pected, the outer portion of the profiles are greatly 
different due to the higher wake velocity through the 
rotor disk in the 2 blade case. 

Consider now a popular algebraic expression used 
to approximate the vortex profile (ref. 28), 

- r » r 

Ve 2tt (r? n +r 2 ")!/" 1 j 

where T* denotes the entire circulation for the vortex. 
To compare the above vortex model with the present 
data, it is assumed that F* has been reached at twice 
the core radius (recall fig. 16). Focusing on data for 
the basic blade at ip = 100°, it appears that the model 
(assuming r* = 2 r c ) is in fairly good agreement when 
n=2 (fig. 18), which is consistent with conclusions 
about this model reported in earlier studies (refs. 28 
and 29). 

Basic Blade Wake- The distribution of vortic- 
ity in the wake of the basic blade provides an imme- 
diate view of the convection path of the vortex dur- 
ing the first revolution of the blade (fig. 19). The 
strength of the vorticity carried by counter-clockwise 
rotating fluid, uj z < 0, is indicated by the level of blue 
saturation in the contour plot (with the trailing vor- 
tex being dominant), while vorticity of the opposite 
sense is represented by various levels of red satura- 
tion. Also included are markers locating the 5 high- 
est neighboring values of vorticity on both extremes. 
These markers serve to expose the irregular complex- 
ion of the vorticity field as well as the migration of iso- 
lated concentrations of vorticity. The pattern shows 
the spiraling accumulation of vorticity from the blade 
wake toward the center of the vortex. The appear- 
ance of positive values of vorticity, u> z > 0, initially 
found along the inboard portion of the wake, are shed 
from the blade where the lift is increasing with blade 
radius. 

The structure of the flow behind the basic blade 
at different wake ages is shown in figs. 20 - 27, and the 
more significant characteristics are summarized in fig. 
28. The location of the center of vorticity appears to 
be quite random from blade revolution to revolution, 
with the extent of the wander increasing with wake 
age. The center of vorticity (based on the centroid of 


high values) is generally at a slightly different location 
than the center of swirl. The midpoint between the 
velocity peaks, based on either v vs. x or u vs. y 
profiles across the vortex, is generally not located at 
the center of swirl. As the vortex ages, the maximum 
values of vorticity and swirl velocity decrease while 
the core radius increases. Tracking the location of the 
trailing vortex based on the centers of swirl at each 
wake age, the rate of movement toward the rotor axis 
of rotation, u, decreases toward an asymptotic value 
while the rate of movement in the direction of the 
flow through the rotor disk, v y increases toward an 
asymptotic value. 

More traditional displays of the vortex structure 
are fashioned in figs. 29 - 36 for the quantities v ) w } 
and c o z in terms of contour plots and cross sections 
through the center of the vortex. The dominant level 
of vorticity is clearly contained in the trailing vor- 
tex, however the peak value of u> z is oftentimes not 
reached monotonically. The out-of-plane component 
of velocity, w, readily exposes the wake deficit, espe- 
cially during the early wake ages (ip < 100°). Maxi- 
mum excursions in w appear to coincide with the peak 
circumferential velocity (as suggested by the profile 
for v) and are directed back toward the rotor blade 
(in keeping with the remainder of the wake deficit). 
Based on the above algebraic model for v$, it has been 
shown (ref. 28) that for n — 2 (which best fits the 
data) that w reaches peak values at about 0.76r c . A 
distinctive zone near the center of the vortex exists 
where u> comes close to matching the velocity in the 
outer 1/r region. 

In certain cases (ip = 10° and 30°) contour plots 
show that w actually points in the opposite direc- 
tion, coinciding with the direction of u> z . It also 
appears that during the early stages of wake devel- 
opment (ip < 100°) that the location of this region 
of “excess” w is offset from the centers of swirl and 
vorticity. Accompanying the display of the v profile 
are portrayals of the velocity relative to the move- 
ment of the vortex in the wake. If the objective is 
to obtain a more symmetric profile in the outer 1/r 
region (thereby rendering a Lagrangian-like appear- 
ance), then the construction based on the convection 
velocity (v — v c , solid line, calculated using actual 
vortex locations) is much more satisfactory than the 
construction based on the velocity measured at the 
interior of the vortex (v — Vi, dashed line). 

0° Vortex Generator- In contrast to the wake 
for the basic blade, a map of the vorticity field with 
the vortex generator attached and fixed at 0° reveals 
the presence of an additional zone of vorticity with 
the same sense as the primary vortex (fig. 37). In 
relation to the primary vortex, this smaller vortex 
persists as a distinct body of fluid for ip < 30° and 
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appears to orbit the primary vortex at an angular 
rate of about 6° for every degree of wake age. Prior 
to the merging of these two vortices, features relating 
to the primary vortex are less certain due to distor- 
tions of the flow field caused by the proximity of the 
smaller vortex (figs. 38 - 45). This is supported by 
a theoretical study of two synthetic vortices which 
shows that the resulting swirling flow becomes very 
different as the vortices begin to overlap (ref. 30). 
A portion of the chordwise vorticity that normally 
would have gone immediately into the primary vor- 
tex may actually have been shed with that from the 
vortex generator. Later wake ages have the appear- 
ance of a single vortex, however they are more intense 
(magnitude of u ) z has increased and the core size has 
decreased) than those measured for the basic blade 
alone (fig. 46). Furthermore, this configuration yields 
an ultimate trailing vortex that has traveled a shorter 
distance (in both x and y) from the tip path traced by 
the basic blade alone, thereby increasing its interac- 
tion with the rotor blade during the next revolution. 
Contour plots and cross sections through the center 
of swirl for the quantities v ) w } and u> 2 are shown in 
figs. 47 - 54. When compared with the basic blade 
case, the u ) z profile shows that the primary vortex is 
initially weaker during the early stages of its devel- 
opment (ip < 30°), but quickly increases in strength 
after the two vortices have merged (figs. 47 - 54). 
The w component of velocity appears to be generally 
weaker than for the basic blade, however, the shape 
of the profile is similar. 

-f-5°,+10 o ,+15° Vortex Generator- The wake 
effects that result from swinging the trailing edge of 
the vortex generator toward the tip of the rotor blade 
are shown in figs. 55 - 108. These angles tend to 
produce tip vortices that have the same sign as that 
trailing from the tip of the rotor blade. As the lift 
on the vortex generator is increased, the tip vortex 
that it produces appears to orbit the primary trailing 
vortex more rapidly than in the 0° case and finally 
merge with the primary vortex at an earlier wake age 
(figs. 55, 73, 91). In each case, distortions in the 
v component of velocity occur during the early wake 
ages through ip — 30° (figs. 56 - 58, 74 - 76, 92 - 94). 
A noticeable jump in the peak magnitude of the w z 
profiles occurs at ip ~ 60°, which coincides with the 
merging of these two vortices (both rotating in the 
same direction). Once the wake has reached ip = 60°, 
all three configurations yield a single vortex with a 
higher peak u) z and a smaller r c than occurred for 
the basic blade alone (figs. 64, 82, 100). Placing the 
vortex generator at +5° and +10° causes the ulti- 
mate trailing vortex to travel farther away from the 
tip path plane (increased y) and with less contrac- 
tion (decreased x) than for the basic blade case. The 


advantage at +5° and +10° of an increased separa- 
tion distance from the rotor plane is not achieved at 
-f 15° possibly due to excessive flow separation from 
the generator at this high angle. 

—5°, — 10°,— 15° Vortex Generator- The wake 
effects that result from swinging the trailing edge of 
the vortex generator away from the tip of the rotor 
blade are shown in figs. 109 - 160. These angles tend 
to produce tip vortices that have the opposite sign as 
that trailing from the tip of the rotor blade. Vorticity 
contours indicate that a secondary vortex does not 
clearly appear until the vortex generator is rotated to 
— 15° (figs. 109, 127, 145). This would suggest that 
the local angle of incidence is different from the geo- 
metric angle of the generator, which is to be expected 
in light of the skewed direction of the approaching 
flow moving around the tip of the rotor blade. Al- 
though not obvious in the contour plots, there is a 
gradually increasing influence on the primary vortex 
as the generator angle is increased. This influence is 
evident in the comparatively flat appearance (which 
increases with the generator angle) on the left side of 
the v profiles from ip ~ 2° to ip = 30° (figs. 110 - 112, 
128 - 130, 146 - 148). Although the sense of the sec- 
ondary vortex is opposite that of the primary trailing 
vortex, it orbits the primary vortex in the same di- 
rection as in the positive vortex generator cases, but 
at half the angular rate. The wake eventually adjusts 
to the presence of the secondary vortex by forming a 
single vortex with a higher peak u> z and a smaller r c 
than occurred for the basic blade alone (figs. 118, 136, 
153). Placing the vortex generator at —5° causes the 
ultimate trailing vortex to travel farther away from 
the tip path plane (increased y) and with less con- 
traction (decreased x) than for the basic blade case. 
The vortex in the —15° case also travels farther away 
from the tip path plane, but now the contraction is 
increased as well. This combination of movements 
places the vortex even farther away from the trace of 
the rotor tip. 

Free Vortex Generator- Assuming that the 
center of gravity is located on the pitch axis, the an- 
gle of this element is free to respond only to local 
flow conditions. Recalling the previous results with 
and without a vortex generator, a review of the con- 
tour maps and cross sectional plots suggests that this 
case is unique (figs. 161 - 178). There is no evidence 
of a secondary vortex present in the contour maps 
(fig. 161), nor are there any distortions in the v pro- 
files during the early stages of wake development (figs. 
162 - 164) comparable to those observed when the 
generator was placed at fixed angles. However, the 
initial vorticity field defining the trailing vortex does 
have an annular appearance similar to the negative 
generator cases, but remains less focused for a longer 
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period of time (fig. 161). This may explain why the 
magnitude of the peak vorticity and the maximum 
swirl velocity are noticeably reduced during the early 
wake ages (fig. 170). The trailing vortex in this case 
has moved a shorter distance from the tip path of the 
rotor as compared to the basic blade case. 

Turbulence Generator- Although this element 
was oriented normal to the blade surface and intended 
to function much like a spoiler, the local flow may 
have recognized it as being more like a flat plate at 
high incidence (and therefore a source of lift). Refer- 
ence to the preceding vortex generator cases will help 
to understand the unusual results for this configura- 
tion (figs. 179 - 196). The contour maps of vorticity 
(fig. 179) show that after the trailing vortex has been 
initially distorted by the turbulence generator, it re- 
mains distributed over a much larger spatial region 
than was observed for all other cases. There does ap- 
pear to be a neighboring zone of positive vorticity at 
— 2° and ^ = 10° (similar to that observed in the 
— 15° generator case), and is undoubtedly the result 
of lift on this element. In keeping with the broadened 
footprint of the resulting trailing vortex, the corre- 
sponding peak velocity and vorticity excursions are 
substantially reduced (fig. 188). The wake contrac- 
tion and its convection away from the rotor plane are 
both less than observed for the blade alone case. 

Overview- The merits of each configuration are 
summarized in terms of vortex strength and location. 
Tables 4 (near field) and 5 (far field) indicate how 
these quantities have changed relative to the basic 
blade. The value of the near-field results may be- 
come more important as the spacing between the ro- 
tor blades is reduced. 

A quantity increases compared to basic blade 
V quantity decreases compared to basic blade 
o no significance difference 


Table 4: Near Field (^ = 30°) 


Con fig 

Max — uj z 

Max vq 

r c 

X 

y 

Free VG 

V 

V 

A 

0 

o 

-15° VG 

A 

A 

O 

0 

o 

-10° VG 

A 

A 

V 

o 

o 

-5° VG 

O 

O 

o 

o 

o 

0° VG 

A 

0 

V 

o 

0 

+5° VG 

o 

V 

o 

o 

o 

+ 10° VG 

o 

V 

o 

o 

V 

+ 15° VG 

V 

V 

o 

o 

V 

90° TG 

V 

V 

A 

o 

V 


In terms of reducing the strength of the trailing 
vortex in the near field (^ = 30°), the free vortex 


Table 5: Far Field (ip = 280°) 


Config 

Max —(jj z 

Max Vq 

r c 

X 

y 

Free VG 

A 

V 

V 

o 

V 

-15° VG 

A 

V 

V 

A 

A 

-10° VG 

A 

A 

V 

A 

O 

-5° VG 

A 

A 

V 

V 

A 

0° VG 

A 

A 

V 

O 
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A 

A 

V 

V 

A 
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A 

V 

V 

A 

+15° VG 

A 

A 

V 

V 

V 

90° TG 

V 

V 

A 

V 
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generator and the turbulence generator are clearly 
superior to the other configurations. However, in 
the far field (ip = 280°), only the turbulence genera- 
tor provides a sustained reduction in vortex strength. 
In terms of increasing the displacement of the trail- 
ing vortex from the path of the rotor tip, none of 
the configurations studied satisfies this requirement 
in the near field. On the other hand, a few cases 
suggest some benefit in the far field, either in terms 
of a greater contraction or a greater convection away 
from the rotor plane. The convection velocity, along 
with the size and location of the trailing vortex, are 
summarized for all cases at ip = 280° in figure 197. 
Relative to the location of the trailing vortex from 
the basic blade, vortices from all other configurations 
remain somewhat clustered nearby and with a maxi- 
mum separation of about one core diameter (or about 
one blade thickness). 

CONCLUSIONS 

1. A robust procedure was developed for analyz- 
ing a vortex dominated flow field, extracting its prin- 
cipal features, and performing a conditional ensemble 
average that accounts for vortex wander. Conditional 
averaging is necessary in order to preserve the inten- 
sity and size of the vortex. 

2. In vortex flows, the highest out-of-plane ve- 
locities often occur where the circumferential velocity 
is highest. When the delay time between PIV pulses 
is too long, interpolations will render in-plane peak 
velocities that are too low. As a result, the vortex 
will appear weaker due to a reduced circumferential 
velocity and an enlarged core size. 

3. The torque produced at 870 RPM in all of 
the vortex generator cases where the angle was fixed 
was lower (up to 6%) than for the basic blade alone. 
However, there was an 8% increase in torque in the 
free (unrestrained) vortex generator case and an 18% 
increase in the turbulence generator case. 

4. The peak value of vorticity generally occurs 
close to the center of swirl. For the basic blade, the 
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component of velocity that is aligned with the axis 
of the vortex (the out-of-plane component, in this 
study) has a somewhat annular shape, but is offset 
from the center of swirl. Peak values of w are di- 
rected back toward the rotor blade and are roughly 
associated with the core diameter defined by the peak 
circumferential velocity, v&. In some cases, the flow 
in the center of this annular region may be pointed 
in the opposite direction (away from the rotor). Dur- 
ing the initial 280° of azimuth the diameter of the 
vortex core increases from about 1.2 to 1.5 times the 
thickness of the rotor blade. 

5. A vortex trailing from a vortex generator with 
a fixed angle will orbit the primary tip vortex at an 
angular rate of about 6° for every degree of wake age 
when the two vortices have the same sense. When the 
subordinate vortex has the opposite sense, the direc- 
tion of the orbit remains the same, but the angular 
rate is greatly reduced. In both cases the subordinate 
vortex cannot be distinguished from the primary vor- 
tex after ip — 60°, and eventually, a more intense 
primary vortex is produced. As the angle of the vor- 
tex generator becomes more negative (producing a 
stronger counter-rotating vortex), the annular shape 
of the w component of velocity becomes more pro- 
nounced as the magnitude of the flow in the interior 
increases in the direction of the rotor. 

6. When the vortex generator is unrestrained, 
the strength of the primary vortex during the initial 
wake ages is significantly reduced. Unlike most of 
the vortex generator cases at fixed angles, this case 
shows no evidence of a subordinate trailing vortex. 
Although the advantages in terms of maximum vor- 
ticity and core radius are diminished after reaching 
ip = 280°, the maximum swirl velocity remains lower 
for all wake ages. 

7. The most dramatic changes occurred in the 
turbulence generator case. The absence of an annular 
shape and the significant reduction in vortex inten- 
sity suggest the effectiveness of the turbulent action 
toward homogenizing the flow inside the trailing vor- 
tex and contributing to its rapid diffusion. Although 
the torque increased by 18%, the maximum vorticity 
was reduced by 65%, the maximum swirl velocity was 
reduced by 57%, and the core size almost doubled 
after reaching ip = 280°. 

8. All of the configurations had some effect on the 
final position of the trailing vortex, with the greatest 
departure from the basic blade case being about one 
blade thickness farther away from the trace of the 
rotor tip. 
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Figure 1: Hover Test Chamber with single-bladed rotor at Ames Research Center. 



Vortex Generator Turbulence Generator 


Figure 2: Single-bladed rotor configuration showing location of tip mounted devices. 
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Figure 3: Concepts studied for altering the trailing vortex. 
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Figure 4: Perspective effect of particle displacement with Stereoscopic PIV. 
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Figure 7: In-situ calibration target as viewed by camera 2. 



Figure 8: Example of image pair with calculation grid. 
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Figure 9: Definition of coordinate system showing origin on rotor- tip path. 
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Ensemble average velocity field 
based on vorticity centroids 


Locate center of swirl pattern 


Calculate core radius 
based on spline fit 


Determine vortex convection 
velocity & path 
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Figure 10: Analysis sequence for extracting trailing vortex characteristics. 
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a : possible swirl center 
b : neighboring node 



Figure 13: Angles used in the search for the center of swirl. 
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Figure 14: Rotor torque resulting from vortex generator and turbulence generator. 
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Figure 16: Effect of vortex convection on local circulation. 
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Figure 17: Comparison with LV data obtained in an earlier experiment (ref. 9). 
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Figure 18: Comparison with algebraic models of vortex. 


27 



Case: Basic 




• max values < 0 

■ max values > 0 

• max values < 0 

■ max values > 0 


....... 

__ V,-..- xfhrir 



f *- 2° K 


CT"*^ *« 100° 






Tl>= 10 


150 


‘iMe. 


• max values < 0 


max values > 0 


max values < 0 


max values > 0 


>.• ■ s-. 


■ ^ 


• max values < 0 

■ max values > 0 

• max values < 0 

■ max values > 0 

! ‘«I Sf's# 4 !®* 

j/u\ 60° |. 


28 0 











LULU '9 S/LU ‘ 8 * 


-6753 


Case: Basic ip= 2° 
std dev allowed= 1.5 
account for wander: focus 
max u z threshold= 50 % 

Vorticity center wander from mean 



Image Pair 

O 



— 1 1 1 I I 

20 40 60 80 100 

r, mm 



6753 s' 1 


-9 mm 


*■*/ 



>.« r ; rn m 


Total image r=-5.9 m 2 /s 


o 



Swirl center 



x, mm 



Figure 20: Vortex statistics at ^ = 2° for the basic blade case. 
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Figure 24: Vortex statistics at ^ - 100° for the basic blade case. 
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Figure 28: Summary of vortex development from $ = 2° -* 280° for the basic blade case 


37 








min max 



Figure 30: Velocity and vorticity components at tp = 10° for the basic blade case. 
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Figure 33: Velocity and vorticity components at ^ — 100° for the basic blade case. 
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Figure 34: Velocity and vorticity components at xjj = 150° for the basic blade case 
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Figure 36: Velocity and vorticity components at ip - 280° for the basic blade case 
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Figure 38. Vortex statistics at rjj — 2° for the 0° vortex generator case. 
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Figure 42: Vortex statistics at ip ~ 100° for the 0° vortex generator case. 
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Figure 44: Vortex statistics at ip ■=■ 210° for the 0° vortex generator case. 
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Figure 45: Vortex statistics at ^ ~ 280° for the 0° vortex generator case. 


54 



200 -100 0 100 200 


• Basic 

■ 0° Vortex generator 




o r Max (-) Vorticity 


o 

lO 


Max Swirl Velocity 






Swirl Center Location 


Swirl Convection 



Wake Age, deg Wake Age deg 


Figure 46: Summary of vortex development from V- = 2° -»■ 280° for the 0° vortex generator case 
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Figure 48: Velocity and vorticity components at ^ = 10° for the 0° vortex generator case 
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Figure 50: Velocity and vorticity components at 4 > = 60° for the 0° vortex generator case 
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Figure 51 : Velocity and vorticity components at = 100° f° r the 0° vortex generator case. 
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Figure 53: Velocity and vorticity components at ip - 210° for the 0° vortex generator case 
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Figure 54. Velocity and vorticity components at tfo = 280° for the 0° vortex generator case. 
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Figure 58. Vortex statistics at ^ — 30° for the H-5° vortex generator case. 
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Figure 60: Vortex statistics at ^ — 100° for the +5° vortex generator case. 
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Figure 61: Vortex statistics at ip = 150° for the +5° vortex generator case. 
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Figure 64: Summary of vortex development from V’ = 2° ->• 280° for the +5° vortex generator case 
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Case: +5° Vortex generator if/= 10 




Figure 66: Velocity and vorticity components at 0 = 10° for the -f5° vortex generator case. 
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Figure 68: Velocity and vorticity components at V> = 30° for the +5° vortex generator case. 
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Figure 69: Velocity and vorticity components at ^ = 100° for the -f5° vortex generator case. 
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Figure 70: Velocity and vorticity components at = 150° for the +5° vortex generator case. 
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Figure 72: Velocity and vorticity components at xp = 280° for the -j-5 0 vortex generator case. 
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Figure 74. Vortex statistics at *0 — 2° for the 4-10° vortex generator case. 
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Figure 76: Vortex statistics at ^ — 30° for the -{-10° vortex generator case. 
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Figure 78: Vortex statistics at 'ip = 100° for the +10° vortex generator case. 
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Figure 80: Vortex statistics at ip = 210° for the H-10° vortex generator case. 
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Figure 82: Summary of vortex development from rjj = 2° -> 280° for the +10° vortex generator case 
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Figure 83: Velocity and vorticity components at rp = 2° for the +10° vortex generator case. 
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Figure 84. Velocity and vorticity components at ip — 10° for the -1- 10° vortex generator case. 
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Figure 85: Velocity and vorticity components at ip = 30° for the +10° vortex generator case 
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Figure 86: Velocity and vorticity components at %j> — 60° for the -j-10° vortex generator case. 
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Figure 87: Velocity and vorticity components at i> = 100° for the +10° vortex generator case. 
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Figure 88: Velocity and vorticity components at $ = 150° for the -h 10° vortex generator case. 
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Figure 90: Velocity and vorticity components at ip = 280° for the +10° vortex generator case. 
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Figure 91: Locations of ui z extrema at rp = 2° ->■ 280° for the +15° vortex generator case. 
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Figure 92: Vortex statistics at if’ — 2° for the -|-15 0 vortex generator case. 
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Figure 94: Vortex statistics at = 30° for the +15° vortex generator case. 
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Figure 96: Vortex statistics at tp = 100° for the +15° vortex generator case. 
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Figure 97: Vortex statistics at tp = 150° for the +15° vortex generator case. 
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Figure 98: Vortex statistics at ip = 210° for the +15° vortex generator case. 
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Figure 100: Summary of vortex development from x\> = 2° ->■ 280° for the +15° vortex generator case. 
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Figure 102: Velocity and vorticity components at xp = 10° for the + 15° vortex generator case. 
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Figure 103: Velocity and vorticity components at ^ = 30° for the +15° vortex generator case 
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Figure 104: Velocity and vorticity components at ip - 60° for the +15° vortex generator case. 
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Figure 105: Velocity and vorticity components at ip = 100° for the +15° vortex generator case. 
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Figure 107: Velocity and vorticity components at ip — 210° for the 4-15° vortex generator case. 
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Figure 108: Velocity and vorticity components at ip — 280° for the +15° vortex generator case. 
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Figure 110: Vortex statistics at ip = 2° for the —5° vortex generator case. 
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Figure 112: Vortex statistics at ip — 30° for the —5° vortex generator case. 
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Figure 114: Vortex statistics at = 100° for the —5° vortex generator case. 
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Figure 115: Vortex statistics at ^ = 150° for the —5° vortex generator case. 
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Figure 116: Vortex statistics at tp = 210° for the -5° vortex generator case. 
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Figure 118: Summary of vortex development from rj> — 2° ->• 280° for the -5° vortex generator case. 
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Figure 119: Velocity and vorticity components at ij) = 2° for the —5° vortex generator case. 
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Figure 120: Velocity and vorticity components at 4 > = 10° for the -5° vortex generator case. 
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Figure 122: Velocity and vorticity components at ip = 60° for the -5° vortex generator case. 
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Figure 124: Velocity and vorticity components at \p — 150° for the —5° vortex generator case. 
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Figure 125: Velocity and vorticity components at xj) = 210° for the —5° vortex generator case. 
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Figure 126. Velocity and vorticity components at ip — 280° for the —5° vortex generator case. 


135 


Case:-10° Vortex generator 


• max values <0 ■ max values > 0 



• max values <0 ■ max values > 0 



• max values <0 ■ max values > 0 



• max values <0 ■ max values > 0 




• max values <0 ■ max values >0 • max values <0 ■ max values > 0 



Figure 127: Locations of u> z extrema at ^ = 2° -> 280° for the -10° vortex generator case. 
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Figure 128: Vortex statistics at ip — 2° for the —10° vortex generator case. 
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Figure 130: Vortex statistics at = 30° for the —10° vortex generator case. 
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Figure 132: Vortex statistics at — 100° for the —10° vortex generator case. 
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Figure 134: Vortex statistics at 0 = 210° for the —10° vortex generator case. 
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Figure 135: Vortex statistics at 4 > = 280° for the -10° vortex generator case. 
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Figure 136: Summary of vortex development from 4 > = 2° ->■ 280° for the -10° vortex generator case. 
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Figure 137: Velocity and vorticity components at xjj = 2° for the —10° vortex generator case. 
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Figure 138: Velocity and vorticity components at xp — 10° for the —10° vortex generator case. 
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Figure 140: Velocity and vorticity components at ip = 60° for the -10° vortex generator case. 
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Figure 141: Velocity and vorticity components at ^ = 100° for the —10° vortex generator case. 
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Figure 144: Velocity and vorticity components at xp = 280° for the —10° vortex generator case. 
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Figure 146: Vortex statistics at xp = 2° for the —15° vortex generator case. 
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Figure 151: Vortex statistics at ip = 210° for the -15° vortex generator case. 
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Figure 152: Vortex statistics at xp = 280° for the —15° vortex generator case. 
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Figure 153: Summary of vortex development from ^ = 2° 280° for the -15° vortex generator case. 
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Figure 154: Velocity and vorticity components at \p — 2° for the -15° vortex generator case. 
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Figure 155: Velocity and vorticity components at ^ — 10° for th e —15° vortex generator case. 
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Figure 156: Velocity and vorticity components at — 30° for the —15° vortex generator case. 
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Figure 157: Velocity and vorticity components at rp = 100° for the -15° vortex generator case. 
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Figure 158: Velocity and vorticity components at = 150° for the —15° vortex generator case. 
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Figure 159: Velocity and vorticity components at ^ = 210° for the -15° vortex generator case. 
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Figure 160. Velocity and vorticity components at ip — 280° for the —15° vortex generator case 
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Figure 161: Locations of u> z extrema at ip = 2° -► 280° for the free vortex generator case 
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Figure 162: Vortex statistics at ip = 2° for the free vortex generator case. 
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Figure 163: Vortex statistics at ip 
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Figure 166: Vortex statistics at tp = 100° for the free vortex generator case. 
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Figure 167: Vortex statistics at ip = 150° for the free vortex generator case. 
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Figure 168: Vortex statistics at — 210° for the free vortex generator case. 
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Figure 169: Vortex statistics at ip = 280° for the free vortex generator case. 
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Figure 170: Summary of vortex development from i/> = 2° 280° for the free vortex generator case 
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Figure 171: Velocity and vorticity components at ip = 2° for the free vortex generator case. 
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Figure 172: Velocity and vorticity components at = 10° for the free vortex generator case. 
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Figure 174: Velocity and vorticity components at rp ~ 60° for the free vortex generator case. 
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Figure 175: Velocity and vorticit.y components at ^ = 100° for the free vortex generator case. 
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Figure 176: Velocity and vorticity components at ip = 150° for the free vortex generator case. 
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Figure 177: Velocity and vorticity components at ^ = 210° for the free vortex generator case. 
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Figure 178: Velocity and vorticity components at ip = 280° for the free vortex generator case. 
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Figure 180: Vortex statistics at ^ = 2° for the turbulence generator case. 
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Figure 181: Vortex statistics at xp = 10° for the turbulence generator case. 
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Figure 182: Vortex statistics at ip = 30° for the turbulence generator case. 
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Figure 183: Vortex statistics at r/; 


u, m/s y, rnm 

-30 -20 -10 0 10 20 30 -50 0 50 100 150 200 250 


Total image T=-4.2 m 2 /s 



J I 1 L 


-320 -240 -160 -80 0 80 
x, mm 


Swirl center 
j u= 2.2 m/s 



26.2 mm 

1 ! I I I ! 

-150 -75 0 75 150 225 300 

y, mm 


60° for the turbulence generator case. 




Image Pair 


Total image T~- 3 . 1 m 2 /s 


v 0 max= 10.8 m/s 
Core radius— 29.6 mm 


0 

01 


0 


14 m/s 


Within Core 
u a ve=- 0 . 6 m/s 
v ave= 1.2 m/s 



— « 1 1 i i 

20 40 60 80 100 

r, mm 



*-320 -240 -160 -80 0 8 ( 

x, mm 


Swirl center 
v= 0.3 m/s 



radius 
27.3 mm 

1 1 j i i i 

0 -275 -200 -125 -50 25 100 

x, mm 



-150 -75 0 75 150 225 30 ( 

y, mm 


Figure 184: Vortex statistics at — 100° for the turbulence generator case. 
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Figure 185: Vortex statistics at i> = 150° for the turbulence generator case. 
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Figure 186: Vortex statistics at xl> — 210° for the turbulence generator case. 
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Figure 187: Vortex statistics at ^ = 280° for the turbulence generator case. 
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Figure 188: Summary of vortex development from V> = 2° 280° for the turbulence generator case 
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Figure 189: Velocity and vorticity components at rjj = 2° for the turbulence generator case. 
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Figure 191: Velocity and vorticity components at ^ — 30° for the turbulence generator case. 
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Figure 192: Velocity and vorticity components at \p = 60° for the turbulence generator case 
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Figure 193: Velocity and vorticity components at rp = 100° for the turbulence generator case. 
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Figure 194: Velocity and vorticity components at 1 > = 150° for the turbulence generator case 
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Figure 195: Velocity and vorticity components at ip = 210° for the turbulence generator case 
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Figure 196: Velocity and vorticity components at V = 280° for the turbulence generator case 
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Figure 197: Summary of convection velocity, size, and location of vortex when - 280° 
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